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Abstract 

In this paper a detailed simulation of irradiated pixel sensors was used to investigate the effects of radiation 
damage on charge sharing and position determination. The simulation implements a model of radiation damage 
by including two defect levels with opposite charge states and trapping of charge carriers. We show that charge 
sharing functions extracted from the simulation can be parameterized as a function of the inter-pixel position and 
used to improve the position determination. For sensors irradiated to 3> = 5.9 x 10 14 n eq /cm 2 a position resolution 
below 15 /im can be achieved after calibration. 



1. Introduction 

The CMS experiment, currently under construc- 
tion at the Large Hadron Collider (LHC) will in- 
clude a silicon pixel detector [1] to allow tracking 
in the region closest to the interaction point. The 
detector will be a key component for reconstruct- 
ing interaction vertices and heavy quark decays in 
a particulary harsh environment, characterized by 
a high track multiplicity and heavy irradiation. At 
the full LHC luminosity the innermost layer, with 
a radius of 4.3 cm, will be exposed to a particle 
fluence of 3 x 10 14 n eq /cm 2 /yr. 

In order to evaluate the effects of irradiation and 
optimize the algorithms for the position determi- 
nation a detailed simulation of the pixel sensors 
was implemented. In [2] we have proven that it 
is possible to adequately describe the charge col- 
lection characteristics of heavily irradiated silicon 
pixel detectors in terms of a tuned double junction 
model which produces a doubly peaked electric 
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field profile across the sensor. The modeling is sup- 
ported by the evidence of doubly peaked electric 
fields obtained from beam test measurements pre- 
sented in [3] . The dependence of the modeled trap 
concentrations upon fluence was presented in [4] 
and the temperature dependence of the model was 
discussed in [5] . Charge sharing and hit reconstruc- 
tion after irradiation was investigated in [6] us- 
ing the so called "77-technique" . In this paper we 
present an alternative approach to hit reconstruc- 
tion and calibration based on charge sharing func- 
tions extracted from the sensor simulation. 

This paper is organized as follows: The sensor 
simulation is described in Section 2, in Section 3 
the hit reconstruction and calibration is discussed. 
The conclusions are given in Section 4. 



2. Sensor simulation 

The results presented in this paper rely upon a 
detailed sensor simulation that includes the mod- 
eling of irradiation effects in silicon. The simula- 
tion, PIXELAV [2,7,8], incorporates the following 
elements: an accurate model of charge deposition 
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by primary hadronic tracks (in particular to model 
delta rays); a realistic 3-D intra-pixel electric field 
map; an established model of charge drift physics 
including mobilities, Hall Effect, and 3-D diffusion; 
a simulation of charge trapping and the signal in- 
duced from trapped charge; and a simulation of 
electronic noise, response, and threshold effects. 
The intra-pixel electric field map was generated 
using tcad 9.0 [9] to simultaneously solve Pois- 
son's Equation, the carrier continuity equations, 
and various charge transport models. 

The simulated devices correspond to the base- 
line sensor design for the CMS barrel pixel detec- 
tor. The sensors are "n-in-n" devices, designed to 
collect charge from n + structures implanted into 
n- bulk silicon. The simulated samples were 22x32 
arrays of 100x150 ^m 2 pixels. The substrate was 
285 /im thick, n-doped silicon. The donor concen- 
tration was set to 1.2 x 10 12 cm~ 3 corresponding 
to a depletion voltage of about 75 V for an unir- 
radiated device. The 4 T magnetic field was set as 
in the CMS configuration and the sensor temper- 
ature to -10° C. The simulation did not include 
the "punch-through" structure on the n + implants 
which is used to provide a high resistance connec- 
tion to ground and to provide the possibility of on- 
wafer IV measurements. 

The effect of irradiation was implemented in the 
TCAD simulation by including two defect levels in 
the forbidden silicon bandgap with opposite charge 
states and trapping of charge carriers. The activa- 
tion energies of the donor and acceptor traps were 
set to (E v + 0.48) eV and (E c - 0.525) eV, respec- 
tively, where Ey and Ec are the valence and con- 
duction band energy level, respectively [10]. The 
trap densities and the capture cross sections for 
electrons and holes were obtained by fitting the 
model to beam test data as described in [2,4]. The 
simulated irradiation fluence was $ = 5.9 x 10 14 
n oq /cm 2 and the reverse bias was set to 300 V. 



3. Hit Reconstruction and Position 
Calibration 

The simulation has been used to study how the 
sharing of charge among neighboring pixels is af- 



fected by radiation damage. A description of this 
behavior is needed to adjust the pixel hit recon- 
struction algorithm as the detector ages. In what 
follows we refer to the reconstruction of barrel hits 
in the r — <p plane, where the charge drift is af- 
fected by the Lorentz deflection. Simulated tracks 
are perendicular to the sensor plane along the 
r — <fi direction. Hits in the pixel detectors are re- 
constructed by first searching for clusters of pixels 
with signals above the readout threshold of 2000 
electrons. The signals are then summed along the 
rows and columns of the cluster to produce 1-D 
projected signal profiles. In the pixel barrel, the az- 
imuthal Lorentz drift produces 2-pixel and 3-pixcl 
wide projected clusters in the local y-direction 
(global ^-direction) . The charge fraction is defined 
as / = Ql/(Qf + Ql) where Ql is the signal of 
the last projected pixel in the cluster and Qf is 
the signal of the first projected pixel. The charge 
fraction of a number of simulated hits is plotted 
in Fig 1(a) as a function of the track coordinate y 
at the midplane of the sensor. The y-coordinate is 
plotted from the center of one pixel (y = 0) to the 
center of the neighboring pixel (y = 100 /im). The 
two-pixel hits are plotted as open diamonds and 
the three-pixel hits are plotted as crosses. These 
functions y = F 2 (f), F 3 (f) arc quite linear and 
have non-negligible width due to fluctuations in 
the deposited charge. After irradiation, the intro- 
duction of trapping states dramatically increases 
the leakage current and produces space charge in 
the detector bulk. The sign of space charge varies 
from negative near the n+ implant to positive near 
the p+ implant resulting in a doubly-peaked elec- 
tric field [10]. The presence of the doubly-peaked 
field and the trapping of the signal carriers affect 
the charge sharing functions which are shown in 
Fig. 1(b) for a fluence of $ = 5.9 x 10 14 n cq /cm 2 . 
We note that nearly all of the 3-pixel clusters have 
vanished, the shape of the 2-pixel function F 2 (f) 
has become non-linear, and there are now 1-pixel 
clusters near the inter-pixel wall (shown as / = 
points). 

The large radiation-induced change in charge 
sharing that occurs during detector operation re- 
quires the implementation of a calibrate-able hit 
reconstruction algorithm. The algorithm should be 
based upon calibration parameters that can be var- 
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Fig. 1. The azimuthal charge-sharing functions for 2- and 
3-pixel clusters in the CMS pixel barrel for new sensors 
(a) and sensors irradiated to a flucnce of <3? = 5.9 X 10 14 
n cq /cm 2 (b). The fraction of charge found in the last pixel 
as compared with both end pixels is plotted as a function 
of the local coordinate y in microns. 

ied in time to remove biases and to ensure the use 
of a technique that optimizes the use of the avail- 
able information. The charge sharing functions are 
well described by the following expression 

\y c + V2 + (w + yi - y 2 ) ■ f a < / < 1 

v=\ (i) 

[y c + A~y c + ( yi +y 2 )/2 f = 

where: y c is the coordinate of the center of the first 
pixel in the cluster, w is the pixel size, y\ and y 2 
define the / = and / = 1 intercepts, the exponent 
a describes the linearity of the function, and the 
offset A is needed to account for the Lorentz-drift- 



induccd asymmetry when the cluster size is one 
pixel. 

The performance of a reconstruction algorithm 
based upon equation 1 is shown in Table 1 for sam- 
ples of simulated events with new and irradiated 
sensors ($ = 5.9 x 10 n eq /cm 2 ). The means and 
the RMS widths of the differences between the re- 
constructed and true hit positions are tabulated 
as functions of the total cluster charge. The large 
cluster charge events are likely to have energetic 
delta rays that spoil the resolution. It is clearly im- 
portant to develop a parameterization of the reso- 
lution for use in track reconstruction. The optimal 
hit reconstruction parameters used in equation 1 
vary considerably before and after irradiation. In 
particular, the exponent a for the two-pixel clus- 
ters varies from 1.00 before irradiation to 0.575 af- 
ter irradiation. 

The parameters of the pixel hit reconstruction 
algorithm are sensitive to bias voltage, operating 
temperature, and irradiation fluence. One can an- 
ticipate that different regions of the detector will 
require different parameter sets and that all will 
change in time. We believe that it will be possible 
to use data to tune the PIXELAV simulation and 
then use the simulation to derive the parameters 
to track the detector aging effects. However, it is 
also possible to extract the shapes of the charge 
sharing functions F(f) directly from data using the 
traditional "ry-technique" [6] . The technique works 
as follows. It is assumed that a group of pixels is 
uniformly illuminated by a set of parallel tracks. 
The traditional formulation uses the variable r/ = 
(Ql-Qf)/(Qf + Ql) = 2/ -1 and with the uni- 
form illumination assumption, it permits the ex- 
traction of F(f) = F(r]) up to an unknown inte- 
gration constant C, 

dN_ _ dN_dy _ Ndy 

dr\ dy dr] w dr] 

■n 

w f , . dN 

^V-NjV-fi+C. (2) 

-l 

This provides a significant constraint on F(f) but 
still requires input from the simulation or another 
technique to determine C. The constant C is re- 
lated to Lorentz-drift which has traditionally been 
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Cluster Charge 


Fraction 


New Sensor 
~Ay RMS 


Irr. Sensor: 


w/o calibration 
RMS 


Irr. Sensor: with calibration 
~Ay RMS 


Q/Qavg < 0.7 


2% 


0.1/im 5.4/im 


26/im 


14/^m 


O.O^m 


7.7/im 


0.7 < Q/Qavg < 1-0 


62% 


0.1/im 7.7/im 


29^m 


12/im 


0.3/im 


9.4/xm 


1.0 < Q/Qavg < 1.5 


30% 


0.4/im 16//m 


29^m 


17/im 


0.2^tm 


15/im 


1.5 < Q/Qavg 


2% 


19/^m 


34/im 


52/im 


11/im 


53/^m 


All Q 


100% 


1/rni 18/im 


29^m 


19/im 


0.6/im 


18/im 



Table 1 

Simulated bias Ay and resolution (RMS) of the pixel hit reconstruction algorithm for different cluster charge bins before 
and after irradiation to $ = 5.9 X 10 14 n cq /cm 2 . The fractions of the sample in each cluster charge bin are also listed. 



studied by searching for an azimuthal cluster size 
minimum as a function of track angle. Unfortu- 
nately, the small azimuthal angular acceptance of 
a pixel barrel module (~ ±10°) and the large size 
of the Lorcntz angle in the pixels (~ 23°) precludes 
this approach. An alternative method, based on 
the measurement of the Lorentz deflection of long 
clusters was proposed in [3]. This method allows 
the measurement of the Lorcntz deflection and the 
extraction of the electric field as function of the 
sensor depth. 



4. Conclusions 

In this paper a detailed simulation of the silicon 
pixel sensors for the CMS experiment was used to 
study charge sharing and the calibration of the po- 
sition reconstruction after heavy irradiation. The 
simulation shows that position resolution can be 
improved after irradiation by parameterizing the 
dependence of charge sharing on the inter-pixel po- 
sition. The procedure requires the tuning of the 
PIXEL AV simulation along the operation of the 
CMS pixel detector, using data collected with col- 
liding beams. 
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